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Abstract— This paper discusses coupled analysis of SAR 
distributions and temperature distributions using an anatomical 
simulation model reconstructed from 2D medical images. The 
distributions calculated by 3D FEM using a large resonant 
cavity applicator for non-invasive deep hyperthermia treatment 
are presented. This applicator is made for treatment of 
abdominal tumors. 


Some electromagnetic therapeutic devices for hyperthermia 
treatment have been used in clinics. However, these applicators 
come with some disadvantages. The abdominal region is 
covered with fat tissue which is resistant to deep heat 
penetration. Furthermore, electromagnetic energy concentrates 
on the convex surfaces of the human body, such as the chests and 
the buttocks. 


In previous studies, it was found that the proposed heating 
system using the resonant cavity applicator was effective for 
heating brain tumors and other small objects. However, when 
heating the abdomen with the developed applicator, areas such 
as the neck, arm, hip or chest were unintentionally heated. 
Therefore, the resonant cavity applicator to overcome these 
problems has been improved. 


To evaluate the effectiveness of the proposed methods, 
coupled analysis of specific absorption rate (SAR) distributions 
and the changing temperature profiles were conducted. First of 
all, an anatomical 3D calculation model from 2D medical images 
was constructed. 


Second, the proposed heating system was presented. A 
cylindrical shield made of an aluminum alloy was installed 
inside the cavity. It was designed to protect healthy tissue from 
concentrated electromagnetic fields. Next, in order to 
concentrate heating energy on deep tumors inside the human 
body, a water bolus was installed around the body. 


Third, the length of the lower inner electrode was adjusted in 
order to control the heating area. To evaluate the effectiveness of 
the proposed methods, temperature distributions were 
calculated by FEM with the 3D anatomical human body model. 

From these results, it was confirmed that the improved 
heating system was able to non-invasively heat abdominal deep 
tumors. 


Index Terms— SAR, Resonant cavity, Anatomical 
model, Finite Element Method (FEM). 


I. INTRODUCTION 


Hyperthermia treatment is based on the clinical fact that a 
tumor is weaker than healthy tissue at the temperature of 
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Fig. 1. Illustration of the heating system. 


42-43°C and can be eliminated by a series of approximately 
one-hour-long heat treatments. A variety of heating methods 
have already been proposed to heat deep tumors. 


Some examples of heating methods are radio frequency 
(RF) capacitive heating applicators and microwave heating 
applicators. Some of these applicators have been in practical 
use. However, all of these heating methods have some 
disadvantages. An successful heating method has not yet been 
realized [1], [2]. 


Therefore, I proposed a non-invasive heating method as 
shown in Fig.1. In this system, a large resonant cavity with 
inner electrodes was used for heating deep tumors in the 
abdominal region of the human body[3]-[5]. 


The human body consists of many organs that have 
different electrical and thermal properties and also have 
various shapes and sizes. Therefore, it was expected that the 
electromagnetic fields would concentrate on the convex parts 
of the human body, such as the chest and the buttocks etc. 


To overcome these problems, I proposed a new heating 
method of using the resonant cavity applicator with 
cylindrical shields made of an aluminum alloy and a water 
bolus. In this new heating system, the human body is covered 
with the cylindrical shields, except for the area to be heated. 


The human body is placed in the gap between the two inner 
electrodes. The surface of the human body is then covered 
with the water bolus to concentrate the electromagnetic 
energy on the deep-seated tumors. 
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Using FEM, I calculated SAR distributions of an 
anatomical human body model, including the bones, muscle 
tissue, fat tissue, lungs, liver, stomach and other internal 
organs. 


II. METHODS 


A. Reconstructing the 3D anatomical human model 


Figure 2 shows the process of reconstructing the 3D human 
body model. The proposed reconstruction method consists of 
3 steps. 


The first step is to collect MRI images taken at intervals of 
several millimeters. In this study, the 2D images of 320 sheets 
were used. After that, shown in Fig.2(a), I traced the outlines 
of all the human tissue. 


The 2nd step is to collect and combine the outlines of all the 
human tissue and make solid models shown in Fig. 2(b). The 
commercial 3D computer-aided design (CAD) software, 
Rhinoceros® (Robert McNeel & Associates) was used in this 
study. 


The last step is to create the FEM model shown in Fig.2(c). 


I used the commercial pre-processor, Hyper Mesh® (Altair 
Engineering, LTD). 


The dimensions of this anatomical human model are 
shown in Fig. 3. This female FEM model has 14 organs. In 
this study, a tumor (3 x 4 x 3 cm) set in the liver is selected as 
the object to be heated. The dimensions of the model are 165 
cm in height, and 40.5 cm in width at the shoulders. The 
average mesh sizes and the electrical parameter values at 
130MHz for each organ are listed in Table I [6]. 


B. The improved resonant cavity applicator 


Figure 4 shows an improved heating system. In Fig. 4, a 
human body is placed in the center of the inner electrodes and 
is heated with the electromagnetic field patterns which are 
generated inside the cavity. No contact is made between the 
human body and the applicator. The elliptic cylindrical shield 
is connected to the cavity wall. 


This shield was designed to protect non-tumorous areas in 
the human body from the electromagnetic heating energy. In 
order to concentrate heating energy on deep tumors, the water 
bolus is set on the human body surface without contacting the 
applicator. 


Figure 5 shows a cross-sectional view of the cavity. The 
cavity is 140 cm in diameter and 100 cm in height. To 
concentrate the heating energy in the center of the cavity, the 
inner electrodes are used. The upper inner electrode is 10 cm 
in diameter and 25 cm in height. The length “L” of the lower 
inner electrode was changed from 15cm to 35cm. In order to 
accommodate differing lengths of the human body, ellipses 
(55 x 30 cm) are located at the head and leg regions. These 
allow the head and leg to rest outside of the cavity. 
Electromagnetic shields are 60 cm in length and are connected 
to the cavity wall. The water bolus, 45 cm in width, 25 cm in 
height and 16 cm in thickness, covers the surface of the human 
body. The water bolus is filled with distilled water. In order to 
heat the liver tumor the human body and bolus were put in a 
position shifted 5 cm from the center of the cavity. 
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Fig. 2. Process of reconstructing the 3D human body model. 
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Fig. 4. Illustration of the improved resonant cavity. 
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(a) Top view of the improved resonant cavity. 
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(b) Side view of the improved resonant cavity. 


Fig. 5. Cross-sectional view of the improved resonant cavity. 


C. Coupled analysis of SAR and temperature distributions by 
3D FEM 


The SAR distribution inside the human body can be 
calculated by equations (1) to (4): 


V?E+k?E=0 (1) 
k? = weu (2) 
1 2 
W, = 77E] (3) 
1 
SAR = 5 Wh (4) 


where Æ is the electric field vector, œ the radial frequency, € 
the dielectric constant, u the magnetic permeability, W, the 
heating power generated inside a human body, ø the electrical 
conductivity, and p the volume density of tissue. Equations (1) 
and (2) can be solved numerically by the FEM [5]-[10]. The 
electrical parameter values at 130MHz for each organ are 
listed in Table I [11]-[14]. 


Using the result data of SAR distributions inside a human 
body, the tissue temperature (T) can be calculated by 
equations (5) to (7) [9], [10], [15]: 


pc = Vk VT +p-SAR-We (5) 
We = (FP )tissue * (PC) blood * (T — Th) (6) 
where c is the specific heat of tissue, « the thermal 


conductivity, JT the temperature of tissues, t the heating time, 
W. the cooling energy by blood flow, F the blood flow rate of 


TABLE I. Electromagnetic Properties and Elements Sizes of Tissues at 










































































130MHz 

Tissue o p Average mesh 
[S/m] [kg/m3] | size [cm] 
Air 0.0 1.00 1.165 5.0 
Bladder 0.30 21.82 1000 2.0 
Bone 0.18 26.21 1790 3.0 
Colon 0.71 76.28 1000 2.0 
Fat 0.036 5.91 900 2.0 
Heart 0.77 83.89 1000 1.0 
Kidney 0.86 89.14 1000 1.0 
Liver 0.51 63.98 1000 0.5 
Lung 0.0 1.00 1.165 2.0 
Muscle 0.72 63.36 1000 1.0 
Sman 1.70 | 87.50 1000 2.0 
Spleen 0.84 82.46 1000 2.0 
Stomach 0.91 74.73 1000 2.0 
Uterus 0.96 75.19 1000 2.0 
Pancreas 0.80 66.67 1000 2.0 
Bolus 0.0 75.00 1000 1.0 
Tumor 0.6 65.00 900 0.2 
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Fig. 6. FEM model (Human body with the applicator). 


each tissue, and T, the temperature of blood. In this study, Tp 
is fixed as 37°C. The initial temperatures of each tissue set to 
37°C, and that of air is 26°C. Heating time was set to 60min. 
Equations (1) and (5) can be solved numerically by the FEM. 
The thermal properties for each organ are listed in Table I 
[16]-[19]. 

Figure 6 shows a finite element mesh for calculating the 
SAR distribution. This FEM model, created by the proposed 
method, and consisting of nonlinear elements, is included in 
the analysis area (250 x 100 cm). The total number of 
elements is 1,898,908. We carried out FEM analysis with the 
3D model using a personal computer constituted with the Intel 
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TABLE II. Thermal Properties of Tissues. 
























































e c K F (Blood flow rate) 
[J/kg/K] [W/m/K] [ml/min/gm] 
Air 1010 0.025 - 
Bladder 3553 0.53 0.31 
Bone 2700 0.22 0.1 
Colon 3012 0.38 1.24 
Fat 2524 0.24 0.21 
Heart 3720 0.54 312 
Kidney 1046 0.52 4.0 
Liver 1050 0.51 3.17 
Lung 3625 0.44 1.14 
Muscle 3634 0.55 0.027 
Small intestine 3012 0.375 2.09 
Spleen 3603 0.54 3.60 
Stomach 3553 0.53 0.53 
Uterus 3542 0.54 0.31 
Pancreas 3601 0.54 0.55 
Bolus 4200 0.60 - 
Tumor 3437 0.50 0.25 

















core 17-4790K and 32GB memories. Here, JMAG (JSOL Co. 
Ltd, Japan), which is a FEM application, was used in 
computer simulations. 


HI. RESULTS 


Figure 7 shows the results of normalized SAR distributions 
calculated by 3D FEM with two applicators, before and after 
the improvements. In these results, the lengths of both inner 
electrodes were 25cm. Figure 7(a) shows the SAR distribution 
with the cavity applicator before the improvement. Here, the 
normalized SAR is given by, 

we (S Snin) (7) 

(S S 
where Sy is the normalized SAR, Smin 1s the minimum SAR, 
Smax 1s the maximum SAR and S is the variable SAR in the 
human body. The resonant frequency was 123.0 MHz. The 
heating energy was concentrated on the hip and the head 
regions. However, in Fig. 7(b), the SAR distribution was 
concentrated on the selected regions with the proposed 
applicator. The resonant frequency was 131.5 MHz. The 
heating power is only concentrated on the gap between the 
inner electrodes. 


max min ) 


In order to reduce the heating energy concentrated on the 
back side, I adjusted the “L”, length of lower inner electrode, 
from 10 to 30cm in 5cm increments. Figure 9 shows the 
results when L is 25cm and 10cm. From these results, the 
concentration of heating energy in the human's backside is 
reduced by adjusting the lower electrode's length. 
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(a) Cavity before the improvement. (Resonant Frequency: 123.0MHz.) 
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(b) Improved cavity. (Resonant Frequency: 131.5MHz.) 


Fig. 7. SAR distributions of side cutting plane. (when L=250mm) 
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Fig. 8. Close up view of SAR distributions. 
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Fig. 9. Close up view of SAR distributions by changing the lower electrode’ s 
length. (when L=25cm and 10cm) 
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Fig. 10. Temperature distributions (when L=10cm). 


Fig.10 shows the results of temperature distribution when L 
is 10cm. From this result, it was found that the heating energy 
was concentrated on the targeted liver tumor. The tumor was 
heated to a maximum of 45°C. 


IV. DISCUSSIONS 


In order to evaluate the effectiveness of the improvements, 
the normalized SAR distributions of before and after the 
improvements were calculated. From Figs. 7 and 8, using the 
previous applicator, the electromagnetic heating energy 
concentrated on convex regions of the human body surface, 
such as the chest, head, neck and hip. However, after the 
improvements, the electromagnetic heating energy does not 
concentrate on the region protected by the shields. The 
heating energy was deeply concentrated on the liver tumor. 


Figure 11 shows the normalized SAR profiles along X-axis 
when L is 25cm. Before the improvements, the normalized 
SAR value of the tumor region was half of the abdominal 
surface region. After the improvements, the maximum 
normalized SAR value is concentrated on the tumor. From 
these results, it was shown that the maximum heating energy 
was deeply concentrated on the liver tumor. However, the 
heating energy was also concentrated on the back region of 
the human body. 


To overcome this problem, the length of the lower inner 
electrode was adjusted to concentrate heating energy on the 
tumor. From Fig. 9, the electromagnetic energy concentrated 
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Fig. 11. Normalized SAR profiles on the X-axis when L=25cm. 
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Fig. 12. Normalized SAR profiles. 


(When changing the length of lower inner electrode.) 


on the back side was changed by adjusting the length of the 
lower inner electrode. 


Figure 12 shows the normalized SAR profiles when 
changing the length “L” of the lower inner electrode. The 
length of the upper electrode was 25cm constantly. All 
profiles show that the heating energy concentrated on the liver 
tumor the most. When the length of electrode becomes shorter, 
the heating energy of the back region decreases. When L is 
15cm, the heating energy of the back region is reduced by 
approximately 15% compared to when L is 30cm in length. 
From this result, the heating energy in the back region can be 
eased by adjusting the length of the lower electrode. Therefore, 
heating energy is concentrated on the deep tumor. And these 
profiles show that the normalized SAR distributions inside of 
the tumor are almost the same although the length of the inner 
electrode was changed. However the SAR distributions of 
back side of the liver was changed because electromagnetic 
properties of these tissue were almost the same. 
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Fig. 13. Temperature profiles on the X-axis.( when L=10cm) 


From this result, it is considered that the SAR distribution 
inside the tumor is not depend on the size of the liver or length 
of inner electrode. 


Coupled analysis of SAR and temperature distributions was 
carried out. From Fig.10, the maximum temperature inside the 
targeted liver tumor was 45°C. 


The surface of the human body was kept cool by the water 
bolus, which was set at a fixed temperature of 28°C. 


As seen in Fig. 13, the temperature distribution results show 
that the tumor region was heated to a maximum of 45°C, 
whole of which is well over the minimum requirement for 
effective hyperthermia treatment (42-43°C). 


From these calculations, it was shown that the improved 
applicator can be effective in heating targeted regions for 
hyperthermia treatment. 


V. CONCLUSION 


This paper discusses coupled analysis of SAR distributions 
and temperature distributions using an anatomical simulation 
model reconstructed from 2D medical images. And I evaluate 
the improved large resonant cavity applicator for 
non-invasively hyperthermia treatment of deeply seated 
abdominal tumor. The proposed heating system was designed 
to protect healthy tissue from concentrated electromagnetic 
fields. 


From SAR distributions, the electromagnetic energy was 
concentrated on the targeted liver tumor. Furthermore, it was 
found that the heating energy of the backside region can be 
eased by adjusting the length of the lower electrode. 


And from the estimated temperature distribution, the liver 
tumor was able to be heated over 43°C using the improved 
applicator. It was confirmed that the improved heating system 
was able to non-invasively heat abdominal deep tumors. 
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